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A BRIEF  REVIEW  OF  SOME  RECENT  RESULTS  THAT  CAN  IMPROVE 
OUR  THEORETICAL  UNDERSTANDING  OF  MAGNETIC  FIELD  RECONNECTION 
AND  THERM AI.IZATION  AS  APPLICABLE  TO  SOLAR  FLARES 

First,  we  should  like  to  remind  the  reader  of  certain  properties  of  the 
tearing  mode;  then  we  shall  attempt  to  relate  the  tearing  mode  developed 
by  Furth,  Killeen,  and  Rosenbluth  (FKR)  (1963)  with  the  vork  of  Sweet 
(1958a, b),  Parker  (1963)  and  Petschek  (1964). 

The  linear  tearing  mode,  as  originally  developed  by  FKR,  shows  that 
the  tearing  mode  can  occur  in  any  magnetic  topology  which  possesses 
shear.  The  mode  was  found  to  be  well  localized  around  surfaces  sometimes 
called  singular  surfaces,  resonant  surfaces,  resistive  layers,  or  mode 
rational  surfaces  where  the  quantity  k*B  vanishes,  k being  the  wave  vector 
of  the  perturbation  and  B the  equilibrium  magnetic  field.  One  should 
notice  that  a large  number  of  surfaces  can  be  excited  if  a spectrum  of  k 
exists  and  |kj  a<l,  where  a ®=(  VI  B|  / |B| )"  * . One  should  further  note  that 

because  more  than  one  surface  can  exist  these  surfaces  can  interact 
as  will  be  discussed  later.  Also  note  that  only  one  component  of  B 
need  vanish  where  k -B  vanishes , i . e . , the  total  B need  not  vanish. 

As  an  example,  applicable  to  a current-carrying  arch  or  filament  (Spicer, 
1976),  the  quantity  k'B  'mB^/r  - <dnBj/R  = B^/r  (m  - nq(r)). 
changes  sign  at  the  mode  rational  surfaces  q = m/n,  where 
q = ir  Bj/B^  R.  B,*,  and  Bj  are  the  poloidal  and 

the  toroidal  components  of  B,  respectively,  and  R the  global  curvature 
of  the  arch.  At  these  surfaces  the  tearing  mode  forms  m x-type  neutral 
points  and  m magnetic  islands  which  result  in  current  perturbations  and 
current  redistribution  and  the  formation  of  current  maxima  and  minima 
in  the  poloidal  direction.  (Figure  1).  Further,  the  tearing  mode  is  really 
a flute  mode  since  where  k-B  vanishes  the  components  of  the  J x B 
restoring  force  vanish,  thereby  permitting  the  continued  growth  of  the 

•The  use  of  the  term  "mode  rational  surface"  is  usually  used  when  discussing  closed  field  configurations,  such  as 

Tokamaks. 

Note  Manuscript  submitted  February  17,  1977. 


1 


mode.  The  fact  that  the  tearing  mode  is  a flute  mode  is  important  when 
one  looks  for  mechanisms  that  can  excite  the  tearing  mode  (Spicer.  1976) 
and  thereby  trigger  a flare. 

As  is  well  known,  one  of  the  first  applications  of  the  tearing  mode  to  the 
flare  phenomenon  was  by  Sturrock  (1°66),  in  his  now  well-known  helmet 
streamer  model  There  are,  however,  a number  of  theoretical  difficulties 
with  this  model.  Basically,  the  tearing  mode  by  itself,  cannot  release 
enough  magnetic  energy;  only  the  energy  contained  initially  in  the  single 
sheet  is  converted,  as  pointed  out  by  Priest  (1975),  since  no  new  magnetic 
flux  is  inwardly  convected.  Note  that  Sturrock  used  only  one  sheet; 
however,  he  could  have  used  a number  of  sheets  by  invoking  a sheared 
field.  Since  no  flux  is  convected  in,  Sturrock  suggested  that  the  tearing 
mode  developed  non-linearly  into  the  Petschek  mechanism.  Is  this 
a valid  assumption?  Kaw  (1976)  appears  to  have  shown  that  it  is.  To  see 
this  we  must  look  at  recent  non-linear  studies  of  the  tearing  mode. 

Rutherford  (1973),  following  a suggestion  by  Kadomstev  and  Pogutse 
(1970),  showed  that  a quasilinear  modification  of  the  basic  current  profile 
6 JOz  * vyl  Bxl/»?  in  sheet  geometry  (Figure  2)  associated  with  the  growth 
of  the  tearing  perturbation  leads  to  a third  order  non-linear  force 
6 JOz  Bxl/c  which  dominates  the  linear  inertial  force  pdV  iy/ dt 
in  the  momentum  equation.  Rutherford  then  replaced  the  inertial  force 
term  with  the  nonlinear  force,  proceeded  with  the  normal  mode  analysis 
and  found  the  similarity  solution  Bx)  * (k  B^0)  *1^  A' ^ t^, 
where  Byo  = dByo/dx  and  A'  is  the  jump  condition  in  the 

logarithmic  derivative  of  Bxj  at  Js.  • B = 0 ; so  that  non-linearly 

the  tearing  mode  grows  as  t^  rather  than  exponentially.  This  result 
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can  be  derived  by  using  some  simple  arguments  but  first  we  will  contrast 
the  difference  between  the  Sweet- Parker  model  and  Petschek's  model 
for  use  later. 

The  Sweet -Parker  model  shows  that  for  a steady  state  configuration, 
the  magnetic  field  convection  produced  by  the  inward  moving  fluid  must  be 
balanced  by  the  resistive  diffusion;  assuming  incompressibility  of  the 
fluid.  Sweet  and  Parker  found 

u = ( V Va/4  it  L)  1 

Thus  in  the  Sweet- Parker  model  resistivity  dominates  the  rate  at  whiefy 
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the  magnetic  field  can  be  taken  into  the  neutral  regions  and  hence  the 
rate  of  magnetic  field  thermalization. 

Petschek  showed  that  the  field  thermalization  problem  is  more 
properly  treated  as  a field  reconnection  problem.  His  results  are  as 
follows: 

. Field  reconnection  signifies  field  annihilation  also,  because  B and 
plasma  kinetic  energy  suffer  changes  at  the  shock  fronts  envisioned  by 
Petshcek. 

. Any  fluid  flow  up  to  u * .1  Va  is  permitted. 

. Resistivity  is  only  important  in  a small  region  about  the  x-type  point 
whose  size  adjusts  itself  to  accommodate  any  flow  (up  to  u 56  .1  VA) 
provided  by  the  boundary  conditions. 

. Field  thermalization  is  not  limited  by  resistivity  but  by  the  fluid 
flow  velocity. 

Following  Kaw  (1976),  we  align  several  Petschek  boxes  end  to  end 
to  form  a periodic  magnetic  island  structure  where  the  reconnection  rate 
and  therefore  the  growth  of  the  magnetic  island  may  be  constrained  either 
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by  resistivity  or  by  fluid  flow  towards  the  singular  layer.  Using  this 
periodic  island  structure  we  first  relate  the  non-linear  results 
of  the  tearing  mode  due  to  Rutherford  (1973)  with  that  of  the  Sweet  - 
Parker  model  both  of  which  are  resistivity  dominated. 

Following  Kaw  (1976)  assume  a black  box  of  length  L and  width 
6 surrounding  the  x-type  neutral  point  (Fig.  (3)).  Plasma  will  flow 
into  the  box  with  velocity  u across  the  interface  having  length  L,  and 
then  flow  along  the  field  lines  into  each  island  and  cause  it  to  expand. 
Assuming  incompressibility,  (which  is  a good  approximation  since 
those  components  of  B^that  do  not  vanish  make  the  effective  pressure 
high  even  though  the  kinetic  pressure  is  small),  the  rate  of  island 
width  growth  is  given  by 


^ » uL/X  , (1) 

dt 

where  X ( =»l/k)  is  the  length  of  the  island,  W the 

island  width  and  u is  evaluated  at  6/2  such  that 

u * E/(  6By0),  (2) 


where  we  will  assume  along  with  FKR  that  resistivity  is  only  important 
around  _k  • fi  = 0 . Using  the  geometrical  relation 

L/6  w X/W  (3) 

and  Maxwell's  equation 


E / X * 3EScl  /dt 


(4) 


we  obtain,  using  (1)  and  (2) 

W * [Bx i X/B^0]1/2.  (5) 
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which  is  the  standard  result  for  an  island  width  (Spicer,  1976; 
Kadomstev  and  Pogutse,  1970).  Ampere's  Equation  together  with 
Ohm's  Law  across  the  resistive  singular  layer  gives 


Jz  = (V  x B)z  = By()/a  * E/tj  * X3Bxi/dt, 

(6) 

since  V x B ^ 0 within  the  singular  layer. 

Using  the  fundamental  assumption  of  Rutherford  that  the  entire 
island  remains  within  the  resistivity  dominated  layer  which  implies 
6 *=  W , equation  (5),  and  the  fact  that  B^.0  % By0/W  within 
the  resistive  layer  we  find  that 


Byi  « BJfO  r}*  t^/a^X  . 


(7) 

we  reach,  along 


Since  X <*>  1/k  and  A'  * 1/a 
with  Kaw,  Rutherford's  result. 

We  can  obtain  the  linear  growth  rate  of  the  tearing  mode  by 
using  the  fact  that  in  the  linear  approximation  the  island  width  is 
much  less  than  the  resistive  layer  width.  Thus  taking  L ^ X and 
using  (1),  (2)  and  (4)  we  have 


dW  % 
dt 


u * E/(B^06) 


hence 


3BX 1 X 

at  (B^q&Y 


W * Bx  i X /(By0  6 ) 

Using 

X/t?  ~Jnr  <By0/a)  81  <By0  w/a> 

we  find 


4 

BX1 


(8) 


Balancing  the  torque  due  to  linear  J x IJ  forces  against 
linear  inertial  forces  gives  6 (ypt?)1/4  (kBJrO)*/2  , 
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where  y = Bx\/Bx\,  so  (8)  yields  the  usual  linear 

growth  rate  . Thus  we  find  for  linear  theory  W « f>  while  for 

quasi-linear  theory  W *=  6 . The  third  possibility  W » 6 

arises  with  plasma  inflow  toward  the  singular  layer. 

So  far  we  have  seen  that  the  Rutherford  theory  is  analogous 

to  the  Sweet -Parker  theory  in  that  both  are  dominated  by  resistivity 

in  the  reconnection  region.  While  this  is  a valid  result  for  a naturally 

evolving  tearing  mode  the  question  arises  - what  occurs  when 

plasma  is  forced  into  the  region  of  reconnection  with  velocity  u.as  is 

the  case  for  the  Petschek  mechanism?  In  this  case,  reconnection  and 

magnetic  island  growth  occur  because  of  reconnection  across  the 

shock  waves  in  the  ideal  region  in  addition  to  the  small  region  around 

the  resistive  layer.  The  island  width  is  then  determined  by  the  flow 

velocity  u given  by  the  external  boundary  conditions  so  that  the  island 

rate  of  growth  is  again  given  by 

dW/dt  * uL/X  , 

where  L is  the  width  of  the  shock  interface.  Again  assuming  L ^ X 
we  obtain 

(9) 

W «=  ut 

so  using 

W * (Bxi  X /By0) 1 ^ 

we  find  (10) 

Bxl  = By0  u2  t2/X  - 

thus  island  width  is  determined  by  the  external  boundary  conditions  on 
the  flow  velocity  and  not  the  resistivity. 
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In  an  arch  model,  as  proposed  by  Spicer  (1976),  this  flow  velocity 
would  be  in  the  radial  direction  so  that  twisting  of  the  feet  of  the  arch 
or  the  propagation  of  a torsional  wave  along  the  arch  should  result  in 
such  a fluid  flow,  thereby  increasing  the  rate  of  field  thermalization  over 
over  that  discussed  by  Spicer  (1976),  and  without  having  to  invoke  the 
phenomenon  of  resonant  overlap  to  be  discussed  below. 

The  results  presented  above  support  the  conclusion  of  Sturrock  (1966), 
i.e. , the  non-linear  evolution  of  the  tearing  mode  by  the  Petschek 
mechanism,  may  occur. 

Since  the  island  widths  can  exceed  the  resistivity -dominated  width 
if  fluid  dominated  and  more  than  one  singular  layer  can  occur  in  a 
sheared  field,  the  question  naturally  arises  "What  would  occur  if  these 
additional  singular  layers  were  to  interact,  i.e. , overlap  one  another.  ?" 

This  question  can  be  answered  by  recalling  that  the  tearing  mode 
is  an  instability  that  requires  destroying  flux  surfaces;  that  is,  the 
adiabatic  invariant  of  ideal  MHD  theory  is  violated,  or 

i 0>  (ID 

dt 

where  <t>  = f B • dS  . This  can  be  made  clearer 

for  an  arch  model  by  converting  the  field  line  equations  in  cylinderical 
geometry 

dr  _ rd0  _ dz  (12) 

Br  ' B*  " Bz 

into  canonical  form.  Imposing  helical  symmetry  through  the  variable 
{ = kz  - m0  , using  V.  • B = 0 , and  using  the 
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stream  function  0 


we  obtain 


and 


1 b<j/ 


Br  = r d£ 


(13) 


R - 1 

* r dr 


where  d\p  / d r = m - kr  B 


(14) 


dr  ’ 

and  1 1/  represents  the 

magnetic  flux  through  a helical  ribbon  defined  by  the  magnetic  axis 
and  a helix  of  constant  phase  £ and  r . Using  p = r^/2 
and  d 0/dt  = (k  - d£/dz)/m  the  field  line 

equations  become 


dp  _ 3^ 
dt  d£ 


(15) 


and 


d£  _ _ 90 
dt  dp 


(16) 


where  0 respresents  the  Hamiltonian,  p the  conjugate 

coordinate  and  £ the  conjugate  momentum,  and  we  have  used 

dz/dt  = B . The  adiabatic  invariants  are 
z 

J = j,  rZ  d £ /2  (17) 
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and 


d 0 <)Jf  0 (18) 

dt  " w = JJ  ' 

J represents  the  flux  tube  area  and  u>  the  frequency  of 
rotation  of  the  flux  line  around  the  flux  tube. 

We  note  that  (15)  and  (16)  can  be  treated  as  equations  of  a non- 
linear oscillator  with  singular  points  9i///9r  = 0 and 

d\p/d(  i.e.  _k  • B_  = 0 and  Br  = 0 .respectively 

Thus  _k  • B = 0 represents  a spatial  resonant  surface  and  in  the 
case  of  an  arch  k_  • B = 0 occurs  when  the  pitch  of  matches 

the  pitch  of  the  perturbation.  Indeed  w can  be  written  as 


As  is  well  known,  a non-linear  oscillator's  frequency  is  amplitude 
dependent.  Because  of  this,  a non-linear  oscillator's  frequency  varies 
when  a resonant  perturbation  acts  on  it,  due  to  the  frequency  dependence 
of  the  driving  force  amplitude.  Thus  the  amplitude  and  frequency  of  the 
oscillator  will  undergo  beating  over  a finite  width  Aco  and 

Ab,  where  w is  the  frequency  and  b the 

amplitude  of  the  driving  perturbation.  In  the  case  of  the  tearing  mode 
Aco  represents  the  width  of  a magnetic  island  and  Ab  the  spread 
in  the  current  perturbation  magnitude  caused  by  the  reconnection  process. 
If  we  denote  e as  the  perturbation  amplitude  due  to  the  tearing  mode  it 
can  be  shown  that  the  width  of  the  resonance  at  • B = 0 , i.e.  , 

the  island  width,  is  given  by 

W = ( e dji  /dr) 1 ^ , 
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which  is  identical  to  our  earlier  result,  where  ^ /rBz  and 

dp /dr  is  the  shear. 

We  can  now  answer  the  question  - what  happens  when  one  or  more 
resonances  interact  or  overlap  one  another?  A definition  of  resonant 
overlap  is 

S = W/d, 

where  W is  the  width  of  one  of  the  resonances  and  d the  distance 
between  them.  When  S < 1 no  interaction  occurs  (Fig.  ( 4))  and  when 
S > 1 interaction  occurs  (Fig.  ( 5 )).  When  S ;>  1 there  is  a 

substantial  increase  in  reconnection  because  islands  are  being  pushed 
against  one  another.  In  addition,  this  interaction  results  in  some  other 
interesting  phenomena.  Suppose  the  resonances  are  well  separated  so 
that  they  do  not  interact.  In  this  case  we  can  superimpose  the  results 
of  a single  resonance.  However,  if  the  resonances  are  close  together 
so  that  S > 1 and  the  phases  of  the  driving  perturbations  are 

random,  field  lines  will  undergo  Brownian  motion  in  their  position,  so 
that  a turbulent  spectrum  of  MHD  waves  may  be  excited. 

Further,  because 

W =«  (e  /dp/dr)1^ 

and  the  distance  between  resonance  surfaces  d is  inversely  proportional 
to  the  shear  it  is  easy  to  see  that  the  strength  of  the  perturbation  required 
for  overlap  becomes  weaker  with  stronger  shear  (Finn,  1975).  We  should 
therefore  expect  strong  overlap  in  arches  with  strong  shear. 
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Recent  non-linear  numerical  simulation  studies  (Schnack  and  Killeen, 
1976)  have  confirmed  that  the  interaction  of  two  singular  layers  leads  to 
a substantial  increase  in  the  degree  of  reconnection  and  the  rate  of 
reconnection.  The  growth  rates  were  found  to  approach  MHD  rates  when 
two  singular  layers  overlapped. 

If  the  width  of  the  islands  is  determined  by  the  flow  velocity  rather  than 
the  resistivity  the  likelihood  of  overlap  occurring  is  also  greater  since  flow 
dominated  reconnection  results  in  larger  island  widths. 

One  can  draw  a number  of  qualitative  predictions  for  flares  if  overlap 
were  to  occur,  they  are: 

. The  volume  of  primary  flare  energy  release  can  be  reduced  because 
of  the  increase  in  effective  reconnecting  volume. 

. Regions  of  strong  magnetic  shear  will  lead  to  flares  with  faster  rise 
times  due  to  the  increase  in  rate  of  magnetic  field  dissipation. 
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Fl|.  2 — Sweet-Parker  mechanism  of  magnetic  annihilation 
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